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ABSTRACT 

Two-dimensional  finite  element  analysis  is  employed  to  predict  the  mechanical  proper¬ 
ties  of  unidirectional  composite  materials  with  isotropic  constituents.  The  composite  is 
treated  as  a  transversely  isotropic  but  homogeneous  continuum,  and  linear  elastic  behav¬ 
ior  is  assumed  for  the  analysis.  Axisymmetric  and  plane  strain  finite  element  formula¬ 
tions  along  with  the  constitutive  relations  for  transversely  isotropic  materials  are  used  to 
obtain  the  apparent  properties  of  the  composite.  The  results  are  compared  with  results 
from  micromechanics  equations  and  experimental  data. 
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INTRODUCTION 


The  evaluation  of  the  mechauicai  properties  of  fiber-reinforced  composite  materials  his 
traditionally  been  a  subject  of  signs Hcant  practical  importance.  Many  analytical  and  numeri¬ 
cal  studies  have  been  performed  on  the  predictioii  of  composite  mechanical  properties.* 

In  particular,  approximate  closed  form  expressions  for  composite  elastic  constants  were 
given  by  Whitney  and  Riley^  based  upon  a  strain  energy  balance  approach  and  classical 
elasticity  theory,  while  bounds  and  expressions  for  the  effective  elastic  moduli  were  given 
by  Hashin  and  Rosen^  based  upon  a  variational  method.  Furthermore,  a  set  of  simplified 
micromechanics  equations  (SME)  to  predict  the  hydro,  thermal,  and  mechanical  properties 
of  composites  have  been  developed  by  Chamis^  based  upon  a  mechanics  of  materials  ap¬ 
proach.  Numerical  approaches  include  finite  difference  methods  by  Tsai,  et  al.^  and 
Adams  and  Doner^  which  were  used  to  predict  composite  elastic  moduli  values.  Finite  ele¬ 
ment  analys's  methods  have  also  been  utilized  to  describe  the  micromechanical  behavior  of 
composites.  Adams  and  Crane^  employed  a  generalized  plane  strain  finite  element  formula¬ 
tion  in  conjunction  with  a  laminated  plate  point  stress  analysis  to  predict  the  stress-strain 
response  of  composite  laminates.  Zhang  and  Evans^  used  an  axisymmetric  finite  element 
approach  and  an  energy  equivalence  principle  to  determine  the  mechanical  properties  of 
composite  with  anisotropic  constituents.  In  the  present  study,  two-dimensional  finite  ele¬ 
ment  analysis  (FEA)  methods  are  used  to  predict  the  mechanical  properties  of  linear  elas¬ 
tic  composites  with  isotropic  constituents.  The  assumptions  and  restrictions  that  have  to 
be  imposed  under  this  type  of  analysis  are  presented  and  discussed.  Numerical  results  are 
also  presented  and  compared  to  those  obtained  from  the  SME  by  Chamis'*  and  experimen¬ 
tal  results  given  by  Adams.^ 

FORMULATION  OF  TNE  BOUNDARY  VALUE  PROBLEM 

The  composite  system  considered  in  the  analysis  is  shown  in  Figure  1.  It  is  assumed 
that  the  system  can  be  represented  by  an  infinite  and  periodic  square  array  of  unidirection- 
ally  oriented  fibers  with  a  relatively  long  axial  dimension  and  of  equal  radius,  r.  The 
assumption  of  periodicity  allows  for  the  isolation  of  a  repeating  volume  element  cell  out¬ 
lined  in  Figure  1  and  shown  in  Figure  2a  for  analysis.  There  are  no  displacements  al¬ 
lowed  across  the  boundaries  of  the  repeating  ceil  and  displacements  are  lestricted  to  those 
that  cause  the  boundaries  to  displace  parallel  to  the  original  boundaries.  Also,  due  to 
symmetry  about  the  y  and  z  axes  only  one  quadrant  of  the  repeating  cell,  the  unit  cell  in 


1.  CHAMIS,  C.  C,  and  SENDECKYJ,  G.  P.  Cmifuf  <m  Theories  Predkw^  Thermoriastk  Proparia  of  Fibrous  Composites.  3.  of  Composite 
Materials,  v.  2, 1968,  p.  332-358. 

2.  WHITNEY,  J.  M.,  and  RILEY,  M.  B.  EUarie  Properties  of  Fiber  Reinforced  Compasiie  Materiab.  AlAA  Journal,  v.  4,  1966,  p.  1537-1542, 

3.  HASHIN,  "L-,  and  ROSEN,  B.  W.  the  Elasdc  Moduli  of  Fiber- Riinfcrctd  Coerporise  Materials.  3.  of  Applied  Mcdiaoics,  v.  31,  1964,  p.  223-232. 

4.  CHAMIS,  C.  C  Simplified  Composite  Mkromechanics  Etftadoea  for  Hygrai,  Thermal,  and  Mechanical  Properues.  NASA  Report  TM-833^, 
National  Aeronautics  and  Space  Administration,  Washington,  DC  1963. 

5.  TSAI,  S.  W.,  ADAMS,  D.  F.,  and  DONER,  D.  R,  Effect  of  Constituent  Material  Pn^enks  on  the  Strengh  of  Ftber-Ranforced  Composite  Ma- 
terials.  AFML-TR-66-190,  U.S.  Air  Force  MaletUls  Labontory  Technical  Report,  1966. 

6.  ADAMS,  D.  F.,  and  DONER,  D.  R.  Transverse  Normal  Loadins  of  a  Unidirectiottat  Composite.  3.  of  Composite  Materials,  v,  I,  1967, 
p.  152-164. 

7.  ADAMS,  D.  F.,  and  CRANE,  D.  A.  Finite  Element  Mkromechtmicai  Analysis  of  a  Urudirectiortal  Composite  Including  Longimdinal  Shear 
Loading.  Cbmputeis  and  Structures,  v.  18, 1964.  p.  1153-1165. 

8.  ZHANG,  W.  C,  and  EVANS,  K.  E,  Numerical  Prediction  of  die  Mecharucal  Properties  of  Anisoinmic  Composite  Materiab.  Computers  and 
Structures,  v.  29,  1988,  p.  413-422. 

9.  ADAMS,  D.  F.  Test  Methods  for  Composte  Materiab.  Seminar  Notes,  U.  S.  Army  Materiab  Technology  Laboratory,  1988. 
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Figure  2b,  need  be  considered  in  the  analysis  to  describe  the  behavior  of  the  repeating 
cell  and,  thus,  completely  characterize  the  state  of  stress  and  strain  of  the  entire  contin¬ 
uum.  Two-dimensional  finite  element  analysis  is  the  method  employed  in  this  study  to 
numerically  solve  the  problem  under  consideration.  Two  separate  hnite  element  mod¬ 
els,  each  employing  two-dimensional  plane  elements,  are  created  and  used  in  the  analy¬ 
sis.  Under  axial  loading  conditions  axial  symmetry  is  assumed  for  the  composite  and  an 
axisymmetric  problem  is  formulated  in  the  i-z  plane.  When  transverse  loading  of  the  com¬ 
posite  is  considered,  plane  strain  conditions  are  assumed,  and  a  plane  strain  mode!  in  the 
2-3  plane  is  employed  in  the  analysis.  Typical  axisymmetric  and  plane  strain  finite  element 
models  are  shown  in  Figures  3a  and  3b,  res^^  jctive*y. 


Figure  1 .  Square  packing  array  unktireciionid  composite. 
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The  appropriate  boundary  conditions  along  the  external  boundarios  of  the  ceil,  as  shown 
in  Figures  3a  and  3b  r^pectively  for  each  problem,  can  be  expressed  as  follows: 


Axisymmetric: 
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Plane  strain: 
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at 

y  =  l,  Z€[0.l]; 
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Plane 

strain  condition; 

u,  =  0 

(10) 

where  ui,  U2,  and  U3  are  the  displacement  components  in  the  1,  2,  and  3  directions,  as  de¬ 
fined  in  Figure  2,  di  and  D2  are  the  naynitudes  of  the  prescribed  displacements  and  d2  and 
D3  are  the  magnitudes  of  the  unknown  displacements.  Furthermore,  continuity  of  tractions 
and  displacements  is  assumed  across  the  fii^r  matrix  interface. 

It  is  assumed  that  the  composite  is  subjected  to  uniform  normal  loads  applied  at  a  dis¬ 
tance  from  the  volume  element  employed  in  the  analysis;  e.g.,  the  remote  stress  field 
shown  in  Figure  1.  When  such  loads  are  considered,  a  complex  nonuniform  state  of  stress  is 
induced  at  the  boundary  of  the  unit  cell  due  to  dissimilar  material  propertimi  of  the  fiber  and 
matrix  phases.  However,  because  of  symmetry,  the  boundaries  of  the  unit  cell  must  displace 
uniformly  under  unifotm  normal  loads  away  from  the  cell  boundaries.  Therefore,  specific  uni¬ 
form  boundary  displacements  in  the  appropriate  direction  can  be  used  to  simulate  the  loading 
conditions  for  the  problem.  Sue!.-  an  approach  applies  the  correct  boundary  conditions,  makes 
the  application  of  loads  easier,  and  accounts  for  the  nonuniforaa  stress  state  at  the  boundary. 
For  the  present  analysis  both  models  are  loaded  using  uniform  prescribed  boundary  dispiace- 
ments,  as  shown  in  Figures  3a  and  3b. 
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ir&c  trerw^isii  irsectioii^cos^ponsats  lat  toB 

celt  wij^e  {H^cfii>ed  dispfau^sMsits  ere  apptied  and  tiie  Earfiefai  ev^afe  strMS  W^ 

can  tie  obtained  front  equillbfittai  cossidenittofis  in  the  appro{»riate  direction  and  can  be  «•> 
press  as: 


I  7*  = 


(11) 


where  S  is  the  cross-sectional  area  of  the  extemai  boundary. 

AastmifiliORa  of  Anaiytis 

In  the  present  tutalysis  an  S2-gia$s/epQiy  system  a  considered.  The  elmdU:  coottant  values 
for  each  distinct  phase  of  the  composite  are  Ihted  in  Table  1.^  Some  of  the  modeling  and 
analysis  assumptions  are  as  follows: 


•  Isotropic  matrix  and  fiber 

•  Composite  is  transversely  isotropic 

•  Linear  clastic  behavior 

•  Continuous  reinforcement  with  perfectly  circular  fibers  of  equal  radius 

•  No  voitfr  in  either  the  fiber  or  the  matrix 

•  Perfect  fiber/matrix  interface 


Table  1.  7VPK>LS2<QLASSFeER  AMD  EPOXY 
MATRIX  PROPERTY  VAUJE8 


Praparty 

S2-Gias8 

Epoocy 

Tansilo  Modulus  (MSI) 

1ZQ0 

a$2 

Poisaon's  Ratio 

0^ 

0.34 

Shaar  Modulus  (MSI) 

5.00 

a23 

In  addition,  effective  stiffness  properties  of  the  composite  are  defined  as  an  average 
measure  of  the  stiffn^  of  the  material  taking  into  account  the  properties  of  all  phases  of 
the  heterogeneous  media  and  their  interactions.  Based  upon  this  definition,  an  averaging 
procedure  is  employed  whereby  the  reaction  forces,  at  the  appropriate  boundaries  of  the  unit 
cell,  along  with  Equation  11  are  used  to  predict  efifective  properties  of  the  composite.  The 
finite  element  analysis  is  used  to  determine  the  displacement  fields  and  reaction  forces  for 
the  composite. 
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Ths  square  packing  array  composite  with  linearly  elastic  components  employed  in  the  anal¬ 
yst  can  be  considered  transversely  botropic.^^  Hie  stress-strain  relations  for  the  composite 
and  its  components  can  be  written  as: 


a,  =Cj  e, 


(12») 


(12b) 


wbei«  Oi  are  the  stress  components,  Sj  are  the  strain  compoi^nts,  and  Qj  and  Sjj  are  the  stiff¬ 
ness  and  compliance  matrices,  respectively.  For  transverse. isotrof^  where  the  2-3  plane  is  the 
plane  of  isotropy  the  normal  components  of  Equation  12b  can  be  written  as; 

£j  =  ^12 ^2  *^12 

£2  “ 'S'j2  tXj +1S22  d^2  "^^23  ^3  (13a) 

£3  —  ^12^1  '^^23^2  '^^^3 

where  in  terms  of  the  elmitic  constants 


and  directions  1,  2,  and  3  refer  to  the  Cartesian  coordinate  system  shown  in  Figure  2. 
When  cylindrical  coordinates  are  considered,  Equation  I3a  can  be  written  as: 


£.  =  S j,<T.  5,2 CT,,  +  5,2 <Ty 

^22  ^  '^23 

£4)  ~5,2<X,  +523  O’,.  ■*’522  00 


(13b) 


where  directions  z,  r,  and  B  refer  to  the  cylindrical  coordinate  system  shown  again  in  Figure  Z 


10.  LEKHNITSKII,  S.  G.  Theay  of  Ekatkity  of  m  Amsotropk  EUuec  Body.  HoMen-Di^,  San  Fnodaoo,  CA,  1963. 
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FiWlfE  ELEM^  IiPl£N^irniTI(»r  - - 

The  finite  element  analysis  method  k  employed  to  solve  the  above  formulated  boundary 
value  problem.  Small  displacements  and  linear  elastic  behavior  are  assumed  throughout  the  so 
lution.  Two  models  each  using  two-dimensional  solid  elements  were  developed  and  u$^  to 
cany  out  the  analysis  and  obtain  all  required  results.  Details  of  these  models  and  how  th^ 
are  used  to  extract  the  necessary  information  are  outlined  below. 

Axlsymmelric  Mod^ 

The  axisymmetric  model  is  used  to  determine  the  behavior  of  the  composite  when 
loaded  in  the  direction  of  the  reinforcement.  The  assumption  of  axial  symmetry  has  been 
shown  to  be  valid  when  considering  a  unidirectional  composite  under  axial  loading;^ ^  and 
the  filament  packing  array  assumed  when  considering  such  a  problem  has  little  effect  on 
the  predicted  property  values.^’ Therefore,  even  though  the  axisymmetric  formulation 
models  a  concentric  cylindrical  geometry  of  the  matrix/fiber  system  as  the  repeating  volume 
element,  the  predicted  property  values  associated  with  axial  loading  of  the  composite 
should  not  be  significantly  affected  by  the  axial  symmetry  assumption.  Furthermore,  the 
computational  and  modeling  advantages  of  axial  symmetry  arc  considerable.  To  determine 
the  displacement  field  and  reaction  forces  for  the  composite  under  axial  loading,  a  uniform 
normal  displacement  is  applied  at  the  2  =  1  boundary  in  the  axial,  2,  direction,  as  shown 
in  Figure  3a.  Note  that  at  the  radial  boundary;  i.e.,  at  r  =  1,  the  radial  stress,  Uf,  van¬ 
ishes  while  the  hoop  stress,  oq,  has  a  finite  value.  Substituting  for  the  elastic  constants 
and  observing  the  conditions  at  the  boundaries  of  the  unit  cell,  the  constitutive  Equations 
of  13b  reduce  to 


Si 


-e 


z 


€2  =e. 


£3 


a,  +- 


(14) 


where  V12  refers  to  the  major  Poisson’s  ratio  of  the  composite  when  loaded  in  the  axial  direc¬ 
tion  and  is  defined  as 


11.  8LOOM,  J.  M.,  wid  WILSON,  H.  B.,  JR.  AxUil  LoaHitg  of  a  Unidevcaonal  Con^poaitc.  J.  of  Compocile  Materiali,  v.  1, p.  268-277. 

' '  ADAMS,  D.  F.,  and  TSAI,  $.  W.  The  Influence  ofRaidom  Filament  Packing  on  the  Dvmveae  Sa/JheB  of  Vnidincaonal  CbmpoMS.  J.  of 
CofflpoMte  Materiib,  v.  3,  1969,  p.  36S-381. 
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(15) 


^12='- 
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Plant  Strain  Modal 

The  plane  strain  model  is  used  to  determine  the  response  of  the  composite  under  trans¬ 
verse  normal  loading.  A  uniform  normal  displacement  is  applied  at  the  =  1  boundary  and 
in  the  transverse,  2,  direction,  as  shown  in  Figure  3b.  Consequently,  a  nonunifonn  displace¬ 
ment  field  is  induced  at  the  2  =  1  free  boundary  of  the  unit  cell,  which  contradicts  bound¬ 
ary  condition  number  7  set  in  the  Assumptions  of  Analysis  Section.  To  e)'iminate  this 
problem  and  verify  boundary  condition  7,  which  traditionally  has  been  imposed  when  similar 
problems  are  sotved,^'^  a  multiceli  approach  is  employed.  The  multicell  approach  utilizes 
multiple  unit  cells  stacked  in  the  z  direction,  while  a  uniform  displacement  load  is  still  main¬ 
tained  at  the  y  =  1  boundary,  as  shown  in  Figure  4.  A  convergence  study  is  performed  to 
determine  the  number  of  cells  needed  for  the  z  boundary  of  the  first  cell  in  the  stack; 
i.e.,  the  cell  boundary'  at  z  1,  to  displace  uniformly  and,  thus,  comply  with  boundary 
condition  7.  Typical  results  of  this  study  are  shown  in  Ingures  S  and  6.  It  is  apparent  from 
these  plots  that  solution  converges  rapidly  and  it  takes  only  four  cells  for  the  boundary 
at  z  »  1  to  displace  uniformly.  Note  that  for  the  multicell  configuration  and  from  equilib¬ 
rium  consiQ»'^i<finns  in  the  z  direction,  the  average  stress  component  o^(z)  =  0  since  there  is 
no  appii  Oh.'  >n  *hat  direction.  Substituting  for  the  elastic  constants  and  observing  the 
conditiors  at  boundaries  of  the  bottom  unit  cell  in  the  multicell  stack,  the  constitutive 
Equations  of  13a  reduce  to  the  following: 


(16) 
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The  numerical  results  obtained  from  the  finite  element  analysis  provide  a  direct 
method  foi  calculating  the  effective  longitudinal  and  transverse  properties  of  the  S2- 
glass/epoxy  composite  system  under  consideration.  The  numerical  procedures  followed  to 
obtain  the  effective  properties  are  reviewed  in  the  Axisymmetrical  Case  and  Plane  Strain 
Case  Sections  below.  Numerical  results  obtained  from  this  analysis  are  compared  to  re¬ 
sults  obtained  from  analytically  developed  micromechanical  equations^  and  the  experimental 
results  of  Reference  9.  All  pertinent  results  obtained  and  comparisons  are  shown  in 
Table  2  and  Figures  7  through  10. 


2  COMPOSTTE  MATERIAL  PROPERTIES 


Method 

%Ftt)er 

VOL(NO 

Ell 

MSI 

E2a 

MSI 

Via 

Vas 

MMA 

0 

0.62 

0.620 

0.340 

0.348 

FEA 

0 

0.62 

0.615 

0.340 

0.337 

MMA 

10 

1.76 

0.88 

0.328 

0.344 

FEA 

10 

1.76 

0.77 

0.325 

0.426 

MMA 

25 

3.46 

1.18 

0.310 

0.340 

FEA 

25 

3.46 

1.04 

0.304 

0.412 

MMA 

SO 

6.31 

1.88 

0.280 

0.330 

FEA 

so 

6.33 

1.92 

0272 

0.306 

MMA 

60 

7.45 

234 

0.268 

0.325 

FEA 

60 

7.44 

2!» 

0.261 

0.255 

EXP 

60 

7.50 

175 

0.280 

— 

MMA 

65 

8.02 

263 

0.262 

0.322 

FEA 

65 

8.01 

3.03 

0.255 

0.^ 

MMA 

100 

120 

120 

0.22 

0.200 

FEA 

100 

120 

11.9 

0.22 

0.218 

NOTE:  MMA  »  Mechanicsl  ol  Materials  Approach;  FEA  »  Fintte  Element  Analysis; 
EXP  s  Experimental  Data. 
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S2 -glass/epoxy  system 
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Figures,  compos  effective  properties. 


AjdsymiiMlric  C«M 

Under  the  condition  of  axial  loading  the  effective  longitudinal  modulus,  and  Poisson’s 
ratio,  Vi2>  the  composite  are  obtained  in  the  following  manner: 

A  uniform  displacement,  di,  is  applied  along  the  z  —  1  boundary  of  the  ax- 
isymmetric  model  and  the  solution  of  the  finite  element  problem  yields  the  un¬ 
known  uniform  displacement,  d2,  along  the  r  1  boundary.  The  average 
strains  in  the  1  and  2  directions  car  then  be  obtained  from 


(17) 


£2 


(18) 


By  employing  Equation  11  the  finite  element  solution  also  yields  the  average  stress  component 
?2;  the  average  hoop  stress  at  the  radial  boundary  is  a  result  of  the  finite  element  solution  as 
well  Using  Equations  14  and  IS  and  substituting  for  the  average  quantities  obtained  from  the  fi¬ 
nite  element  solution,  the  effective  elastic  constants  Vj2  E  j  given  by 


(19) 


(20) 


Plan*  Strain  Caaa 

Under  the  ^ndition  of  transverse  loading,  the  effective  transverse  modulus,  and 
Poisson’s  ratio,  V23  are  obtained  in  the  following  manner: 

A  uniform  displacement,  D2,  is  applied  along  the  y  =  1  boundary  and  the  solu¬ 
tion  of  the  finite  element  problem  yields  the  unknown  uniform  displacement, 

D3,  along  the  2  =  1  boundary.  It  should  be  noted  that  only  the  first  unit  cel! 
in  the  stack  of  the  multicell  problem  is  used  to  obtain  numerical  results  for 
the  plane  strain  case.  The  average  strain  in  the  2  and  3  directions  can  then 
be  obtained  from 
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By  employing  Equation  11  the  finite  element  solution  yields  the  average  stress  component 
9 .  Using  Equation  16  and  substituting  for  the  average  quantities  obtained  from  the  Enite 
element  solution,  the  effective  elastic  constants  ^23  given  by 


(23) 


V  <^2  ) 


(24) 


CONCLUSIONS 

In  this  study,  the  Enite  element  method  was  used  to  predict  the  effective  elastic 
properties  of  an  S2*glass/epo:qr  transversely  isotropic  composite.  Circular  fibers  in  a  square 
packing  array  and  linear  elastic  behavior  for  the  fiber  and  matrix  phases  was  assumed  in 
the  formulation  of  the  boundary  value  problem.  Axisymmetric  and  plane  strain  finite  ele¬ 
ment  models  of  the  problem  under  consideration  were  used  simultaneously  to  btain  the 
elastic  constants  of  the  composite.  Numerical  results  obtained  from  the  finite  element  solu¬ 
tion  were  compared  with  results  obtained  using  classical  micromechanics  equations,^  as  well 
as  with  experimentally  obtained  results  at  60%  Eber  volume.^  It  was  found  that  values  ob¬ 
tained  for  the  longitudinal  modulus  £  ,  were  in  very  good  agreement  with  both  theoretical 
and  experimental  results,  while  values  Tor  the  rest  of  the  elastic  constants  did  not  compare 
very  well.  Numerical  values  of  the  elective  transverse  modulus,  £  ,  were  found  to  be 
greater  than  both  analytical  and  experimental  values,  particularly  at  nigher  Eber  volumes. 
This  can  be  attributed  primarily  to  the  assumptions  made  in  the  analysis  and  the  two-dimen¬ 
sional  formulation  of  the  finite  element  problem.  The  plane  strain  assumption  seems  to  re¬ 
sult  in  a  transversely  stiffer  composite  system.  It  also  produces  erroneous  results  for  the 
effective  in-plane  Poisson’s  ratio,  V23  .  Furthermore,  in  utilizing  two-dimensional  Enite  ele¬ 
ment  analysis  two  problems  had  to  be  solved  to  obtain  all  of  the  constants.  However,  the 
geometry  employed  in  the  two  Enite  element  problems  is  not  the  same.  The  axisymmetric 
problem  solves  the  composite  cylinder  assemblage  problem,  which  has  been  traditionally 


16 


employed  to  obtain  approximate  closed  form  solutions  for  the  composite  elastic  constants, 
while  the  plane  strain  problem  solves  a  single  Hber  system  within  a  square  boundary.  This  dif¬ 
ference  in  geometric  constraints  at  the  boundari^  might  have  affected  the  behavior  of  the 
composite  system,  thus  introducing  errors  in  the  calculation  of  the  elastic  constants.  The 
muiticell  plane  strain  problem  successfully  applied  the  correct  condition  at  the  cell  bound¬ 
ary  at  z  =  1.  Figures  11  through  14  show  how  the  model  behaves  and  bow  the  z  bound¬ 
ary  displaces  as  the  number  of  unit  cells  increases  from  one  to  four. 


17 


RguivU. 


In  conclusion,  it  seems  that  the  approach  employed  to  model  the  micromcchanical  re¬ 
sponse  of  organic  composite  materials  in  this  analysis  is  inadequate.  A  full  three-di¬ 
mensional  finite  element  analysis  study  must  be  undertaken  to  relax  some  of  the 
constraints  and  limitations  imposed  by  the  two-dimecsicnal  fonsuiation.  It  should  also  be 
noted  that  assumptions  made  in  the  analysis,  such  as  the  Hber  packing  array  and  interface 
properties,  must  be  reviewed  and  assessed. 
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